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Abstract 

The price of financial assets are, since Q, considered to be de- 
scribed by a (discrete or continuous) time sequence of random vari- 
ables, i.e a stochastic process. Sharp scaling exponents or unifractal 
behavior of such processes has been reported in several works || || [f| 
H H| . In this letter we investigate the question of scaling transforma- 
tion of price processes by establishing a new connexion between non- 
linear group theoretical methods and multifractal methods developed 
in mathematical physics. Using two sets of financial chronological time 
series, we show that the scaling transformation is a non- linear group 
action on the moments of the price increments. Its linear part has a 
spectral decomposition that puts in evidence a multifractal behavior 
of the price increments. 



1 Introduction 

One of the pillars of modern physics is the covariance of theories under certain 
group actions. What is particular to a given application, such as initial 
and boundary conditions usually breaks the symmetries of the theory. The 
symmetry group of observed data is therefore usually much smaller than 
the covariance group of the theory. An example is hydrodynamics where 
the equations are invariant under space-time translation and scaling, but 
where the solutions are not, in general, invariant. For theories that are 
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covariant under scaling (to be specific, we can think of Navier-Stokes or 
Korteweg-De Vries equation) the situation is clear: the scaling properties, 
such as the spectral decomposition of the solution at each time are given 
by the scaling properties of the initial condition, the boundary conditions 
and external forces. The study of scaling transformation properties in the 
domain of economics and finance is more complicated because the evolution 
equations (or even the theory) governing the dynamics are largely unknown. 
It is thus not possible, in this case, to separate these scaling properties into 
a general property of an underlying theory and into what is particular to 
the situation under study. The observed financial chronological data results 
at least to some extent, of the particularities of each market and not only 
from a general abstract dynamics. Therefore there is no a priori reason to 
expect the data to exhibit simple properties under scaling transformation. 
Keeping this simple observation in mind, we will base our analysis of the 
scaling transformation on methods adapted to physical systems with complex 
behavior: 

i) Multifractal analysis of fully turbulent systems introduced in []7j and within 
that approach further developed inversion techniques developed in || (see 
also |J). 

ii) Non-linear group representation theory developed in JTOj and applied to 
many non-linear evolution equations in mathematical physics (see |TT| for 
recent contributions). 

We apply these methods to two sets of financial chronological series : 

1) Foreign exchange rate DM/$: The data set provided by Olsen and 
Associates contains worldwide 1,472,241 bid-ask quotes for US dollar-German 
mark exchanges rates from 1 October 1992 until 30 September 1993. Tick 
by tick data are irregularly spaced in recorded time. To obtain price values 
at a regular time, we use linear interpolation between the two recorded time 
that immediately precede and follow the regular time. We obtain in this way, 
for a regular time of 15 seconds, 1,059,648 data. Our study focuses on the 
average price which is the mean of the bid and ask price. 

2) Stock index CAC 40: The data set provided by the "Societe de Bourse 
Frangaise" contains 1,045,890 quotes of the CAC 40 index from 3 January 
1993 until 31 December 1996. Tick by tick data are regularly recorded every 
30 seconds, during opening hours (everyday from 10 am until 5 pm except 
week ends and national holidays). Our data base consist of the daily registers 
to which a constant has been subtracted such that the value at 10am is equal 
to the value of the previous day at 5pm. The subtracted jump process (with 
jumps at fixed times) can be analyzed on its own. This separation allows for 
a finer analysis of the rest of the process. 
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Using these data sets, we obtained three new results. First, the scaling 
transformation of the moments of the observed probability distribution is a 
non-linear representation that is well approximated by a linear representa- 
tion, for small scaling parameters. This linear representation turns out to be 
diagonal. Secondly, the function of the order of the moment, defined by the 
spectrum of the generator is (non trivially) concave. This shows, by definition 
@, that the data are multifractal. Note that the concavity in the case of FX 
market (DM/$) can partially be deduced from [12] and is confirmed, inde- 
pendently of our work, by |13| . Our third new result is an explicit expression 
of the family of probability distributions of price increments corresponding 
to different time increments. 

For larger values of the scaling parameter, the linear approximation breaks 
down and the non linear terms of the representation has to be considered. 
The analysis of this letter can also be applied to the SP 500 index, where 
the results should be compared to the (unifractal) scaling behavior found in 
0. It should also be compared with the, from the point of view of finance, 
more fundamental approach of stochastic time transformation (subordinate 



processes) that were applied to SP500 [Q Ul5] . These points are left for 
future investigations. 



2 The mathematical framework 

We suppose that the financial variable is described by a stochastic process 
{u(t)) t >o such that the set of increments u{t + r) — u(t), r > has a well 
defined transformation property under scaling of the time increment r, th 
a r, a > 0. To avoid complications, irrelevant for the quite crude application 
reported in this letter, we suppose that (u(t)) t >o is stationary. Moreover, 
we will only consider the absolute value \u(t + r) — u(t)\ of increments. Let 
w(t) = \u(t) — m(0)|, r > 0. This means that for each (scaling) a > 0, there is 
a map T a of the set W = {w(t)\t > 0} such that T a {w{r)) = w(ar). A group 
action T of the scaling (dilatation) group ID (the set of strictly positive real 
numbers) on the set W is then defined, i.e T a b(x) = T a (Tt>(x)) and T e (x) = x 
for a, b G ID, x G W, where e = 1 is the identity element in ID. In the cases 
under consideration in this letter, it follows from the observed time series 
that the estimated probability distribution p w ( T ) in IR of w(t) is different for 
different r > 0. This is enough to ensure the existence of the action T, and 
moreover shows that T gives a group action T, on the set M. = {p w {t)\t > 0} 
of probability distributions, defined by T a (p w ^) = p w ( aT y 

The group action T is not linear, in spite of its appearance. To ex- 
plicit properties of the scaling action T, we change the coordinates of the 
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elements in A4. As in the case of fully developed turbulence, we use the 
moments as coordinates. For q G Ai, let the moment vector be the sequence 
S(q) = (S r (q)) r >o, where S r (q) = J °° x r q(x)dx and r G IR + . Here we suppose 
that the set M. of probability measures is such that S r (q) exists for all orders 
r and moreover that q is determined by its moments of order r G IR + (which 
is the case if for example the Fourier Transform of elements in Ai are quasi 
analytic). Let S be the image (in the space C(JR + ) of continuous real func- 
tions on IR + ) of Ai under the coordinate transformation S. The image U of 
the group action f is given by U a = SoTaoS' 1 , i.e U a (S(p w ( T ))) = S(p w{ar )). 

In the case U is a linear diagonal representation, it has the form U a = Ua~\ 
where for given real numbers ( r with r G 1R + : 

UP(m) = (a c 'm r ) r > (1) 

for a G ID and m G C(IR + ), m r corresponding to a moment of order r. We 
note that {Cr| r £ 1R + } is the spectrum of the generator of the representa- 
tion U w . When U is a non linear perturbation of there are algorithms 
permitting its construction. However they are outside the scope of this let- 
ter |L0| [p| . For commodity we denote s r (r) = S r (p w f T \) which is the r-th 
component of U T (S(p w m)). An accurate and explicit approximation of the 
inverse transformation S" -1 , of the moment vectors s r (r) to probability distri- 



bution p w r T ) has been developed in p] P |]16| . This permits to obtain directly 
from experimental data, an explicit formula for the family M. = {p w ( t )}t>o 
of probabilities. In fact, for each r G IR + we can determine an element p w r T ) 
by the formulas: 



xp w ( T )(x) = p(lnx) 

<i In s r (r) 

a(r,r) = —^U. (2) 

.... , . . c?lns r (r) 1 . . 1 d 2 \ns r (T) . . 

\np(a(r,r)) = In a r (r) - r £ 2 ^ ~ 2 dr 2 - (3) 



where r G 1R H 



3 Results 

When the representation {7 is linear it follows from expression (0) that 
lns r (r) = A r + Crhir, where A r and ( r are independent of r. Figure ([TJA) 
shows that, in the case of FX DM/$, this is satisfied, to a good approximation 
with time increments r and moments of order r in the interval 11 < r < 2896 
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Figure 1: (1A): lns r (r) vs lnr for r = 1 . . . 10 in the case of FX DM/$ 
index. The scaling law lns. r (r) = A r + £ r lnr displayed as straight lines for 
11 < t < 2896 minutes (delimited by the orange arrow) gives ( r as the slope. 
(IB): ( r determined by scaling law in the case of FX DM/$ index. 




Figure 2: (2A): lns r (r) vs lnr for r = 1 ... 10 in the case of CAC40 index. 
The scaling law lns r (r) = A r + C r In r displayed as straight lines for 1 < r < 
2048 minutes (delimited by the orange arrow) gives ( r as the slope. (2B): ( r 
determined by scaling law in the case of CAC40 index. 



5 



minutes and 1 < r < 10. In contrast, for CAC40 the domain of validity of 
the linear approximation also contains the small values of r: 1 < r < 2048 
minutes and 1 < r < 10 (see figure (0A)). Outside this domain in the (r,r) 
plane, the linear representation approximation breaks down. Inside the do- 
main of validity of the linear representation approximation, the spectrum of 
the generator is presented in figure (QB) (resp. (@B)) in the case of FX DM/$ 
(resp. CAC40). The function r i— > ( r is in both cases (non trivially) con- 
cave, which by definition (see [[!]]) shows that the system has a multifractal 
behavior. 

Finally, we present in figure (|3|A) and figure (QB) (resp. figure (QA) and 
figure (f|B)) probability densities (M.A.M) given by (0) and (|]), for r = 8 
minutes and r = 512 minutes in the case of FX DM/$ index (resp. CAC40). 
In all the cases, the experimental probability distribution is well approxi- 
mated, for a large range of price increments, by the corresponding probabil- 
ity distributions in the family {p w ( t )}t>o constructed by the inverse method 
developed in || and Jy| . Other commonly used probability distributions are 
also presented in the figures for illustration. 
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Figure 3: (3A): Presentation of the probability density function at r = 8 
minutes for FX DM/$ index and comparison with empirical data. (3B): 
Same presentation at r = 512 minutes for FX DM/$ index. 
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Figure 4: (4A): Presentation of the probability density function at r = 8 
minutes for CAC40 index and comparison with empirical data. (4B): Same 
presentation at r = 512 minutes for CAC40 index. 
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